functionFolder = cd('..');
oldFolder = cd('input\');
load('DLCIsoquantWorkspace.mat')
cd('..');

if optfig.plotfig == 1
    
    [fontname,dimension,lw,folder,color,style,marker,markersize,fontsize_tit,fontsize_ax,fontsize_lab,fontsize_leg] = fn_optfig(optfig);
    
    [X,Y] = meshgrid(basisPointImpactGrid',dollarVolumes');
    %surf(X,Y,CCBEuro','FaceAlpha',0.5)
    fig = figure;
    
    % Create axes
    axes1 = axes('Parent',fig);
    hold(axes1,'on');

    % Create surf
    surf(X,Y,TwinShareDivergence','FaceAlpha',0.5);
    view(axes1,[149.699997425033 9.60000030710107]);
    xlabel('$$\alpha$$ Coefficient','interpreter','latex','FontSize',fontsize_lab,'FontName',fontname)
    ylabel('Average Dollar Trading Volume','interpreter','latex','FontSize',fontsize_lab,'FontName',fontname)
    zlabel('DLC Ratio','interpreter','latex','FontSize',fontsize_lab,'FontName',fontname)
    title('\$10B Welfare Gain DLC Isoquants', 'interpreter','latex','FontSize',fontsize_tit,'FontName',fontname)
    ax = gca;
    ax.FontName = fontname;
end

functionFolder = cd('.');
oldFolder = cd('output/');
saveas(gcf,'05_SurfacePlot','epsc')
cd(functionFolder);